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ABSTRACT 

We present the low-resolution spectroscopic and UBVRI broad-band photometric in- 
vestigations of the Type lib supernova 2011fu, discovered in UGC 01626. The photo- 
metric follow-up of this event has been initiated a few days after the explosion and 
covers a period of about 175 days. The early-phase light curve shows a rise followed 
by steep decay in all bands and shares properties very similar to that seen in case of 
SN 1993J, with a possible detection of the adiabatic cooling phase. Modelling of the 
quasi-bolometric light curve suggests that the progenitor had an extended (~ lx 10 13 
cm), low-mass (~ 0.1 Mq) H-rich envelope on top of a dense, compact (~ 2 x 10 11 
cm), more massive (~ 1.1 Mq) He-rich core. The nickel mass synthesized during the 
explosion was found to be ~ 0.21 Mq, slightly larger than seen in case of other Type 
lib SNe. The spectral modelling performed with SYNOW suggests that the early-phase 
line velocities for H and Fen features were ~ 16000 km s _1 and ~ 14000 km s _1 , 
respectively. Then the velocities declined up to day +40 and became nearly constant 
at later epochs. 
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1 INTRODUCTION 

It is commonly recognized that core-collapse supernovae 
(CCSNe) represent the final stag e s of the life of massive stars 
(M > 8-10 MfTi) (|Heger et all 120031 ; lAnderson fc James! . 
2009; ISmarttl,l2009l ). Generally, the fate of massive stars is 
;overn e d by i t s mass, meta l licity, r otation and magne t ic fiel d 
Frverl, 1 19991 ; iHeger et all 120031 ; IWooslev fc Jankal. l2005h . 
Massive stars show a wide variety in these fundamental 
parameters, causing diverse observational properties among 
various types of CCSNe. The presence of dominant hydro- 
gen lines in the spectra of Type II SNe strongly suggests 
that their progenitors belong to massive stars which are still 
surrounded by significantly thick hydrogen envelope before 
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the explosion (see iFilippenkol . 1 19971 . for a review on differ- 
ent types of SNe). On the contrary, Type lb events are H- 
deficient but He is still present in their spectra, unlike to 
Type Ic SNe, where both H and He features are absent. 
After the disc overy of SN 1987K, another class , term ed as 
Type lib (see IFilippenkol . 1 19881 ; IWooslev et all . Il987l h was 
included in the CCSN zoo, and the observational properties 
of these SNe closely resemble those of Type II SNe during 
the early phases, while they are more similar to Type Ib/c 
events at later epochs. 

However, in a few cases, the spectral classification of 
Type lib SNe i s more controversial : for example, SN 2 000H 



jBenetti et al.1 . |2000| ; iBranch et al.1. |2002|; lElmhamdi et al 
200rJ); SN 2 00 3bg l|Filippenko fc Chornock[ [20031 



this 



Soderberg et~all . l2006t ); SN 2007Y (althoug h 
event is classified as Type Ib/c by e.g. iMonar d| l|2007f) ; 
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IStritzinger et alj (|2009l ) however, iMaurer et al.l 
suggested that it is a Type lib) and SN_2009mg (jPrieto 



sugge sted that it is a Type lib) and alN2009mg (|rnetc , 
l2009l : lRoming et alJ . [2009al : IStritzinge3 . l2010l )). SNe of Type 
lib are further divided into two subgroups: Type cllb with 
compact progenitors like SNe 1996cb, 2001ig and 2008ax, 
and Type e llb with extended progenitor s, e.g SNe 1993J 
and 2001gd l|Chevalier fc Soderberg| , l2010l ). 

Type lib and Type Ib/c SNe are collective l y kno wn 
as "stripped envelope" CCSNe |Clocchiatti et all 1 19971 ) as 
the outer envelopes of hydrogen and/or helium of their pro- 
genitors are partially or completely removed before the ex- 
plosion. The possible physi cal mechan i sms b ehind this pro- 
cess may be stellar winds l|Puls et al.l . 12008) or interaction 
with a companion star in a binary syst em where mass trans- 
fer oc curs due to Roche lobe overflow l|Podsiadlowski et all 
Il992l ). There are several studies about the discovery of 
the progenitors of Type lib SNe but the debate about 
how they man age to keep on l y a t hin layer of hydro- 
gen, is still on (lAldering et all KL993: iMaund et all [2004 : 



Ryder et al.l.l2006l;ICrockett et alll2008l:ISonbas et allhoOS 



20111 : IVan Dvk et al 



Arcavi ct al., 2011; Maund ct al 
201 it ISoderberg et all |2012| ). 

To date, approximately 7l{3 Type lib SNe are 
known, but only a few of them have been prop- 
well-st udied. Amo ng them, SNe 
19881 ): 1993 J (ISchmidt et all 



erly mon itored and 
1987K llFilippenkol. 
Il993l : ILewis et alf 
1996cb ilQiu et all 



2009 ]: iMazzali et all 



2008 



Roming et al 




IRichmond et all 119*941) 



2003bg (iHamuv et ai 



2 008ax (iPastorello et al 



Taube nbcrgcr ct al.. 120111 
2011ei (Milisavlicv ic et al 
2011dh jArcavi et aT ~ 



2011 



2009bl : IChomock et all 

2009m g |Oates et al.1 . |2012|) : 
2012|) more recentl 



2011 



Marti- VidaletaL, 201 



Maund et all l201ll: IVan Dyk et all l201ll: iBersten et aU , 
20121: iBTetenholz et all I2OI2I: iHoresh et all 



Krauss et al.1 . 120121 : ISoderberg et all 120121 : IVinko et~al 



20121 ) were remarkably well-studied 
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et all 



An interesting property of the observed light curves 
(LCs) of a few Type lib SNe is the initial peak and rapid 
decline followed by a subsequent rise and a secondary max- 
imum. The first peak is thought to be due to the break- 
out of the SN shock f rom the extended progenitor envelope 
jFalk fc Arnettl 1 19771 ). The properties of the shock break- 
out peak depend on the envelope mass and the density struc- 
ture of the outer layers. The shock break-out phase can last 
from seconds to days. Therefore, early discovery and rapid- 
cadence early-time observations may help in understanding 
the propert ies of the outer envel ope of massive stars in a 
better way l|Gal-Yam et al.l . l201ll ). 

In this paper, we present the results from photomet- 
ric and spectroscopic monitoring of SN 201 lfu starting 
shortly after the discovery and extending up to nebular 
phases. The photometric and spectroscopic properties of 
this event reveal that SN 2011fu is a Type lib supernova. 
The type determination for this SN was verified with SNID 
ijBlondin fc Tonrvl . l2007f ), highlighting an excellent resem- 
blance of the object with SN 1993J. 

The paper is organized as follows. The photometric and 




Figure 1. The \Z-band image of the SN 201 lfu field around the 
galaxy UGC 01626, observed on 2011 November 16 with the 1-m 
ST, India. The SN is marked with a black arrow. The reference 
standard stars used for calibration are marked with numbers 1-8. 
On this image, north is up and east is to the left. 



spectroscopic observations are presented in Section 2, where 
the methods for data reduction and analysis are described. 
In Section 3, we analyze the light and colour curves. In Sec- 
tions 4, we describe and model the bolometric LC. Section 5 
deals with the spectroscopic modelling using the SYNDW code. 
In Section 6, we discuss the metallicites of the host galaxy 
of SN 201 lfu along with those of other CCSNe. Finally, the 
results are summarized in Section 7. 



2 OBSERVATIONS AND DATA ANALYSIS 

SN 2011fu was discovered in a spiral arm of the galaxy UGC 
01626 (type SAB(rs)c) by F. Ciabattari and E. Mazzoni 
jCiabattari et alll201ll ) on 2011 September 21.04 (UT) with 
a 0.5-m Newtonian telescope in the course of the Italian Su- 
pernovae Search Project. The brightness of the SN at the 
time of discovery was reported to be at mag ~ 16 (unfil- 
tered). It was located 2" west and 26" north of the cen- 
ter of the host galaxy, with coordi nates a = 02 h 08 m 21°41 , 
S = +41°29'12'.'3 (equinox 2000.0) (|Ciabattari et all l201ll l. 
The host galaxy has a heliocentric velocity and redshift of 
5543 ± 11 km s" 1 and z = 0.018489 ± 0.00003/0, respec- 
tively. The first spectrum of SN2011fu was obtained on 
2011 September 23.84 UT with the Ekar-Copernico 1.82- 
m telescope (range 3 60-810 nm; resolution 2.2 nm) by 
iTomasella et al.l (|201ll ). showing a blue continuum with su- 
perimposed weak H and He I 587.6-nm features, which led 
to the classification as a young Type II SN. 



2.1 Optical Photometry 

The prompt photometric follow-up of SN 2011fu started 
shortly after the discovery and continued using three 



1 http:heasarc.gsfc.nasa.gov/W3Browse/star- 
catalog/ asiagosn.html 
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ground-based telescopes in India. The majority of the obser- 
vations were made using the 2-m Himalayan Chandra Tele- 
scope (HCT) of the Indian Astronomical Observatory, Hanle 
and the 1-m Sampurnanand Telescope (ST) at the Aryab- 
hatta Research Institute of observational sciencES (ARIES), 
Nainital, India. All observations were performed in Bessell 
UBVRI bands. 

The HCT photometric observations started on 2011 
September 28 using the Himalaya Faint Object Spectro- 
graph Camera (HFOSC). The central 2k x 2k region of a 
2k x 4k SITe CCD chip was used for imaging which pro- 
vided an image scale of 0.296 arcsec pixel -1 across a 10 x 
10 arcmin 2 field-of-view. 

Further photometric observations were carried out using 
a 2k x 2k CCD camera at the f/13 Cassegrain focus of the 1- 
m ST telescope situated at ARIES, Nainital. The CCD chip 
has square pixels of 24 x 24 urn, a scale of 0.38 arcsec per 
pixel and the entire chip covers a field of 13 x 13 arcmin 2 
on the sky. The gain and readout noise of the CCD camera 
are 10 electron per ADU and 5.3 electrons, respectively. A 
finding chart showing the field of the SN 2011fu along with 
the local standard stars is presented in Fig. [T] 

In addition, we also observed this SN in V and R bands 
on 2011 Decemb er 01 and 2012 March 02 using the 1.3m tele- 
scope DFOTEI (|Sagar et all l201ll . I2012T) . recently installed 
at Devasthal, Naintial. 

To improve the signal-to-noise ratio (S/N), all the pho- 
tometric observations were carried out with 2x2 binning. 
Along with science frames several bias and twilight flat 
frames were also collected. Alignment and determination of 
mean FWHM on all science frames were performed after 
the usual bias subtraction, flat fielding and cosmic-ray re- 
moval. The standard tasks available in irafQ and DAOPHOllfl 
(|Stetsorl fl987l . 1 1992h were used for pre-processing and pho- 
tometry. 

The pre-processing steps for images taken with all three 
telescopes were performed in a similar fashion. The stellar 
FWHM on the V^-band frames typically varied from 2" to 
4", with a median value of around 2'.'5. We also co-added 
individual frames, wherever necessary, before computing the 
final photometry. 

For phot ometric calibra tion, we observed the standard 
field PG0231 (|Landoltl 12009) in UBVRI bands with the 1-m 
ST on 2011 December 17 under good photometric conditions 
(transparent sky, seeing FWHM in V ~ 2"). The profile fit- 
ting technique was applied for the photometry of SN 201 lfu 



3 DFOT uses 2048 x 2048 ANDOR CCD camera having 13.5 x 
13.5 urn pixels mounted at the f/4 Cassegrain focus of the tele- 
scope. With 0.54 arcsec per pixel plate scale, the entire chip covers 
a 18 X 18 arcmin 2 field-of-view on the sky. The CCD can be read 
out with 31, 62, 500 and 1000 kHz speed, with system RMS noise 
of 2.5, 4.1, 6.5, 7 electrons and gain of 0.7, 1.4, 2, 2 electron/ADU 
respectively. We selected the 500 kHz readout frequency during 
our observations. 

4 IRAF stands for Image Reduction and Analysis Facility dis- 
tributed by the National Optical Astronomy Observatories which 
is operated by the Association of Universities for research in As- 
tronomy, Inc. under co-operative agreement with the National 
Science Foundation. 

5 DAOPHOT stands for Dominion Astrophysical Observatory Pho- 
tometry. 



and Landolt field and then instrumental magnitudes were 
converted into standard system follow ing least - squar e lin- 
ear regression procedures outlined in IStetsor] |1992). At- 
mospheric extinction values (0.57, 0.28, 0.17, 0.11 and 0.07 
mag per unit airmass for U, B, V, R and I bands, respec - 
tively) for the site were adopted from iKumar et all (2000). 
The chosen Landolt stars for calibration had the bright- 
ness range of 12.77 ^ V ^ 16.11 mag and colour range 
of —0.33 ^ B — V ^ 1.45 mag. Using these stars, transfor- 
mation to the standard system was derived by applying the 
following zero-points and colour coefficients: 

u - U = (7.27 ± 0.01) + (-0.08 ± 0.02)([/ - B) 

b-B = (4.90 ± 0.004) + (-0.04 ± 0.01)(B - V) 

v — V = (4.34 ± 0.01) + (-0.04 ± 0.01)(B - V) 

r-R = (4.19 ± 0.01) + (-0.04 ± 0.01)(V - R) 

i-I= (4.60 ± 0.02) + (0.04 ± 0.02)(V - I) 

Here U, B, V, R, I are the catalogue magnitudes and 
u, b, v, r, i are the corresponding instrumental magnitudes. 
Table [1] lists the coordinates and magnitudes of the eight 
local secondary standard stars in the SN field. 

To estimate the possible contribution from the host 
galaxy to the measured supernova fluxes, we used the ISI^ 
image subtraction package. We acquired deep images (hav- 
ing total exposure times of more than 20 minutes) in BVRI 
bands with the HCT telescope on 25 August 2012 under 
good sky conditions. As the supernova was not detected in 
anyone of these frames, we used them as template frames for 
image subtraction. We found minor differences, not exceed- 
ing 0.1 mag, between the SN magnitudes with and without 
applying the image subtraction for the data at later epochs 
i.e. 70 days after the first observation. The final results of our 
SN photometry (without applying image subtraction correc- 
tions) along with robustly determined PSF errors, are pre- 
sented in Table [21 

2.2 Spectroscopy 

Spectroscopic observations of SN 201 lfu were obtained at 
8 epochs, between 2011 September 28 (JD 2455833.27) and 
December 22 (JD 2455918.11). A journal of these observa- 
tions is given in Table[3] The SN spectra were taken with the 
HFOSC instrument mounted at the 2-m Himalayan Chan- 
dra Telescope. All spectra were obtained using grisms Gr#7 
(wavelength range 3500 - 7800 A) and Gr#8 (wavelength 
range 5200 - 9200 A). FeAr and FeNe arc lamp spectra 
were applied for wavelength calibration. Spectrophotometric 
standard were also observed with a broader slit to correct 
for the instrumental response and flux calibration. 

The reduction of the spectroscopic data were carried out 
in a standard manner using various tasks available within 
IRAF. First, all images were bias-subtracted and flat fielded. 
Then, one dimensional spectra were extracted from the two- 
dimensiona l clean e d im ages using the optimal extraction 
algorithm l|Hornei I1986T I. The wavelength calibration was 
computed using the arc spectra mentioned above. The ac- 
curacy of the wavelength calibration was checked using the 

6 http:/ /www2. iap.fr/users/alard/package. html 
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Table 1. Identification number (ID), coordinates (a, 8) and calibrated magnitudes of standard stars in the field of SN 2011fu. 



Star 


Q J2000 


<5j2000 


U 


B 


V 


R 


/ 


ID 


(h m s) 


r") 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


1 


02 08 19.66 


+41 30 53.4 


15.59+0.02 


15.47+0.02 


14.80+0.01 


14.38+0.02 


14.03+0.02 


2 


02 08 14.25 


+41 27 51.8 


18.39+0.09 


18.17+0.02 


17.59+0.02 


17.16+0.02 


16.83+0.03 


3 


02 08 27.70 


+41 29 11.4 


16.29+0.02 


15.70+0.02 


14.86+0.02 


14.34+0.02 


13.93+0.03 


1 


02 08 26.92 


+41 30 07.5 


17.53+0.04 


17.58+0.02 


17.07+0.02 


16.70+0.02 


16.40+0.03 


5 


02 08 21.94 


+41 29 26.9 


16.78+0.03 


16.77+0.02 


16.21+0.01 


15.84+0.02 


15.54+0.02 


6 


02 08 12.72 


+41 32 34.8 


16.36+0.02 


16.45+0.02 


15.85+0.01 


15.44+0.02 


15.07+0.02 


7 


02 08 08.04 


+41 30 41.3 


16.32+0.02 


16.07+0.02 


15.38+0.01 


14.97+0.02 


14.61+0.02 


8 


02 08 06.54 


+41 30 43.7 


17.23+0.03 


16.03+0.02 


14.93+0.02 


14.32+0.02 


13.78+0.02 
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Figure 2. Observed UBVRI light curves of SN 2011fu. For clarity, the light curves in different bands have been shifted vertically by the 
values indicated in the legend. Black solid lines represent the light curves of SN 1993J over-plotted with appropriate shifts. The explosion 
date of SN 2011fu was taken to be 2011 September 18 ± 2, as described in subsection 3.1. 
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Table 2. Photometric observational log of SN 2011fu 



JD 


Phase a 


U 


B 


V 


R 




Telescope 




(Days) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 





n p" p" o o o no 

2455833.23 


+10.73 


-1 *7 'if' i i \ no 

17.36+ 0.03 


2455834.49 


+11.99 


-1 t / • / 1 i n no 

17.66± 0.03 


2455836.15 


+13.65 




n .1 p" p" o o t n/' 1 

2455837.26 


+14.76 


-1 t to 1 n no 

17.73± 0.03 


o .1 p" p" o .1 -i no 

2455841.23 


+18.73 


-1 t c r" 1 n n o 

17.65± 0.03 


2455842.30 


+19.80 


i T ff\ l n n r" 

17.69± 0.05 


n ^ p" p" o J o nT 

2455843.27 


+20.77 


-1 t s ' o I n nn 

17.62+ 0.09 


2455844.21 


+21.71 


~t ^7 a a \ n n 

17.48± 0.06 


2455845.44 


+22.94 


-it ^ o 1 n n o 

17.43± 0.03 


2455845.43 


+22.93 




2455846.44 


+23.94 


1 t r" -i I n n o 

17.51± 0.03 


2455846.43 


+23.93 




O A P" P" O J n J"l 

2455849.41 


+26.91 


1 t o I n no 

17.68± 0.03 


2455850.40 


+27.89 


17.76+ 0.04 


2455851.31 


+28.81 


i o oi i n -1 o 

18.21+ 0.13 


n .1 p" p" o p" t ori 

2455857.30 


+34.80 


-I o nn 1 n n/ -1 

18.99+ 0.06 


n .1 p" p" o p" o o o 

2455858.32 


+35.82 


-I n nn 1 n nn 

19.00+ 0.09 


2455859.18 


+36.68 




2455860.30 


+37.80 


1 n 1 n I n n o 

19.10+ 0.08 


2455862.30 


+39.80 


19.24+ 0.15 


n .1 p" p" o *i ■! n 

2455864.40 


+41.90 


-1 n p" r" 1 n n o 

19.55+ 0.08 


n .1 p" p" o p" or* 

2455865.35 


+42.85 




2455866.23 


+43.73 




2455866.26 


+43.76 




2455875.22 


+52.72 




2455879.28 


+56.78 




n .1 p" p" o o -i n /"> 

2455881.26 


+58.76 




2455882.33 


+59.83 




2455884.27 


+61.78 


1 n n /i 1 n n p" 

19.94+ 0.05 


2455894.23 


+71.73 


19.67+ 0.20 


o p" p" o n no 

2455896.28 


+73.78 


-I n tp - 1 n nT 

19.75+ 0.07 


2455897.08 


+74.58 




n p" p" o n o on 

2455898.30 


+75.80 




2455900.17 


+77.66 




2455901.24 


+78.74 




2455904.28 


+81.78 




2455909.18 


+86.68 




n p" p" n -i n no 

2455912.28 


+89.78 




2455910. 24 


+90.73 


1 n pt n I n 1 o 

19.56+ 0.12 


n ^ rrm o -i o 

2455918.18 


+95.68 




2455919.11 


+96.61 




o j rrnrio -it 

2455922.17 


+99.67 




n p" p" n n *i -in 

2455924.10 


+101.60 




2455929.17 


+106.67 




2455930.27 


+107.76 




2455930.27 


i ~i n o to 

+108.73 




n p" p" n o n -i p" 

2455932.15 


+109.65 




n p" p" n o z' 1 or 

2455936.25 


+113.75 




2455937.16 


+114.66 




2455938.08 


+115.58 




2455939.23 


+116.73 




2455947.06 


+124.56 




2455947.16 


+124.66 




2455953.08 


+130.58 




2455954.16 


+131.66 




2455963.14 


+140.64 




2455967.11 


+144.61 




2455969.09 


+146.59 




2455976.07 


+153.56 




2455979.12 


+156.62 




2455989.08 


+166.58 




2455998.08 


+175.58 




2456159.34 


+336.84 
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16 
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17. 
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.69+ 





.02 


16. 


19.05 


± 


0.03 


17. 


77+ 





.01 


17 


19.12 


± 


0.02 


17. 


81+ 





.01 


17 


18.97 


± 


0.05 


17. 


.88+ 





.02 


17. 


19.10 


± 


0.02 


17. 


86+ 





.02 


17. 








18. 


06+ 





.06 


17. 


19.30 


± 


0.10 


18. 


17+ 





.03 


17 








18. 


19+ 





.02 


17 


19.38 


± 


0.07 


18. 


.19+ 





.03 


17. 


19.49 


± 


0.03 


18. 


28+ 





.01 


17. 


19.37 


± 


0.05 


18. 


41+ 





.03 


17. 


19.49 


± 


0.03 


18. 


43+ 





.01 


17. 








18. 


39+ 





.03 


17. 


19.28 


± 


0.05 


18. 


41+ 





.02 


17. 


19.47 


± 


0.08 


18. 


50+ 





.03 


17 


19.36 


± 


0.06 


18. 


48+ 





.04 


17. 


19.31 


± 


0.17 










17. 


19.52 


± 


0.07 


18. 


60+ 





.03 


17. 


19.45 


± 


0.05 


18. 


64+ 





.03 


18. 


19.47 


± 


0.05 


18. 


59+ 





.03 


18 


19.55 


± 


0.03 










18. 


19.58 


± 


0.04 


18. 


75+ 





.01 


18. 


19.48 


± 


0.06 


18. 


74+ 





.03 


18 


19.65 


± 


0.03 


18. 


.86+ 





.01 


18 


19.73 


± 


0.09 


18. 


92+ 





,04 


18. 








19. 


04+ 





.11 


18. 








18. 


89+ 





.15 


18. 


19.80 0.16 


19. 


.00+ 





.05 




19.69 ± 0.09 


19. 


09+ 





.04 


18. 


19.70 ± 0.04 


19. 


11+ 





.02 


18. 








19. 


00+ 





.07 










19. 


21+ 





.02 


18. 
















18. 


19.64 ± 0.08 










18. 








19. 


33+ 





.02 


18. 








19. 


23+ 





.07 










19. 


18+ 





.07 


18. 
















18. 








19. 


40+ 





.03 


18. 
















18. 








19. 


.59+ 





.09 


19. 
















19. 



>22.5 >22 



99 


± 





02 


16 


.74 


± 





.02 


HCT 


11 


± 





.02 


16 


.86 


± 





.02 


HCT 


08 


± 





.02 


16. 


88 


± 





.02 


HCT 


02 


± 





.02 


16. 


83 


± 





.02 


HCT 


77 


± 





.02 


16. 


59 


± 





.02 


ST 


70 


± 





.02 


16. 


56 


± 





.03 


b 1 


67 


± 





.02 


16. 


.57 


± 





.04 


b 1 


62 


± 





.02 


16. 


48 


± 





.03 


bl 


59 


± 





.02 


16. 


45 


± 





.03 


ST 


54 


± 





.02 


16. 


43 


± 





.02 


HCT 


57 


± 





.02 


16. 


43 


± 





.03 


ST 


54 


± 





.02 


16. 


45 


± 





.02 


HCT 


53 


± 





.02 


16. 


38 


± 





.02 


HCT 


57 


± 





.02 


16. 


43 


± 





.03 


HCT 


56 


± 





.02 


16. 


44 


± 





.03 


ST 


80 


± 





.02 


16. 


58 


± 





.03 


b 1 


84 


± 





.02 


16. 


61 


± 





.03 


ST 


92 


± 





.02 


16. 


69 


± 





.02 


HCT 


93 


± 





.02 


16. 


68 


± 





.02 


ST 


97 


± 





.02 


16. 


71 


± 





.03 


ST 


09 


± 





.02 


16. 


84 


± 





.02 


HCT 


11 


± 





.02 


16. 


85 


± 





.02 


HCT 


14 


± 





.02 


16. 


83 


± 





.03 


ST 


17 


± 





.02 


16. 


88 


± 





.02 


HCT 


33 


± 





.04 


16. 


95 


± 





.05 


bl 


55 


± 





.03 


17. 


18 


± 





.03 


ST 


50 


± 





.02 


17. 


18 


± 





.03 


HCT 


56 


± 





.02 


17. 


20 


± 





.03 


ST 


62 


± 





.01 


17. 


30 


± 





.02 


HCT 


76 


± 





.02 


17. 


42 


± 





.03 


ST 


87 


± 





.02 


17. 


57 


± 





.03 


HCT 


79 


± 





.03 












DFOT 


79 


± 





.02 


17. 


45 


± 





.03 


b 1 


83 


± 





.03 


17. 


53 


± 





.03 


ST 


88 


± 





.03 


17. 


39 


± 





.03 


b 1 


97 


± 





.05 


17. 


.72 


± 





.04 


HCT 


96 


± 





.03 


17. 


63 


± 





.03 


b 1 


06 


± 





.03 


17. 


67 


± 





.04 


ST 


09 


± 





.03 


17. 


70 


± 





.03 


ST 


21 


± 





.03 


17. 


82 


± 





.03 


HCT 


22 


± 





.02 


17. 


87 


± 





.02 


HCT 


24 


± 





.03 


17. 


81 


± 





.04 


ST 


34 


± 





.02 


18. 


02 


± 





.02 


HCT 


36 


± 





.05 


18. 


00 


± 





.04 


b 1 


47 


± 





.07 


17. 


99 


± 





.12 


bl 


29 


± 





.10 


17. 


93 


± 





.10 


ST 










18. 


11 


± 





.03 


HCT 










18. 


27 


± 





.06 


HCT 


58 


± 





.03 


18. 


13 


± 





.04 


HCT 


65 


± 





.03 


18. 


35 


± 





.03 


HCT 










18. 


11 


± 





.05 


ST 


74 


± 





.02 


18. 


46 


± 





.04 


HCT 


55 


± 





.04 


18. 


16 


± 





.07 


ST 


68 


± 





.03 


18. 


38 


± 





.06 


ST 


87 


± 





.02 


18. 


49 


± 





.03 


HCT 










18. 


54 


± 





10 


ST 


87 


± 





.05 


18. 


.48 


± 





.08 


ST 


91 


± 





.05 












ST 


88 


± 





.03 


18. 


58 


± 





.07 


ST 


92 


± 





.05 


18. 


.51 


± 





.09 


ST 


24 


± 





.10 


18. 


.87 


± 





.08 


DFOT, S 


07 


± 





.05 


18. 


89 


± 





.10 


ST 


>21.5 






>21 




HCT 



6 Brajesh Kumar et al. 

Table 3. Log of spectroscopic observations of SN 2011fu. 



Date 


J.D. 


Phase™ 


Range 


Resolution 






(Days) 


(A) 


(A) 


2011-09-28 


2455833.28 


+10.78 


3500-7800; 5200-9250 


7 


2011-09-29 


2455834.41 


+11.92 


3500-7800; 5200-9250 


7 


2011-10-01 


2455836.23 


+ 13.73 


3500-7800; 5200-9250 


7 


2011-10-14 


2455849.43 


+26.93 


3500-7800; 5200-9250 


7 


2011-10-29 


2455864.41 


+41.91 


3500-7800; 5200-9250 


7 


2011-10-31 


2455866.36 


+43.86 


3500-7800; 5200-9250 


7 


2011-11-23 


2455889.15 


+66.65 


3500-7800; 5200-9250 


7 


2011-12-22 


2455918.11 


+95.61 


3500-7800 


7 



with reference to the explosion epoch JD 2455822.5 



night sky emission lines and small shifts were applied to the 
observed spectra whenever required. The instrumental re- 
sponse curves were determined using the spectrophotomet- 
ric standards observed on the same night as the SN, and the 
SN spectra were calibrated to a relative flux scale. When the 
spectrophotometric standards could not be observed, the re- 
sponse curve based on observations in a night close in time 
to the SN observation was adopted. The flux calibrated spec- 
tra in the two regions were combined to a weighted mean to 
obtain the final spectrum on a relative flux scale. 

Finally, the spectra were brought to an absolute flux 
scale using zero points determined from the calibrated, 
broad-band UBVRI magnitudes. The SN spectra were also 
corrected for the redshift of the host galaxy (z = 0.018), and 
de-reddened assuming a total reddening of E(B — V)= 0.22 
mag (see Sub-section 3.3). The telluric lines have not been 
removed from the spectra. 



3 LIGHT CURVES OF SN 2011FU 

In this section, we present the multi-band light curves of 
SN 201 lfu and their comparison with SN 1993J light curves 
and their temporal properties. A brief discussion about the 
explosion epoch of SN 201 lfu is presented in the following 
sub-section. 



(i) For CCSNe of Type lb and lib, the explosion dates 
have bee n estimated to be ~ 20 d ays prior to t he l^ -band 
maxima (Rich ardson et al. U2006I ; iDrout et all l201ll ) (see 
also iMilisavlievic et al. , 2012). 

(ii) Type lib SNe also exhibit bluer B - V c olour ~ 40 
days after the explosion (|Pastorello et all 120081 ). giving an 
indication about th e explosion epoch . 

(iii) The SNID l|Blondin fc TonrvL l2007h fitting on ini- 
tial four spectra of SN 2011fu indicates that explosion of 
this event would have occurred around 2011 September 20. 
However, SNID fit for the later three epochs of spectra (af- 
ter V band maximum) gives rise to 2011 September 17 as 
the explosion date. 

(iv) In some of the well studied type lib SNe, the explo- 
sion epoch is better constrained (e.g SN 1993J, SN 2008ax 
and SN 2011dh) and their early light curve features indi- 
cate that the adiabatic cooling phase may be observable for 
several days after the explosion and this du ration depends 



upon the volume of the photospheric shell dRoming et al 



2009b) 

10.0.9- 1 



as determ i ned in case of SN 19 9 3J (|Wheeler et al_ 
Il993l; (Lewis et all 11994 iBarbon et ail ll995Tl, SN 2008a x 



l|Roming et alll2009bl ) and SN 2011dh ifArcavi et all 1201 ll ) 

Based on above evidences, we have adopted 2011 
September 18 ± 2 as the explosion epoch for SN 2011fu and 
it will be used for the further discussions in this article. 



3.1 Explosion epoch of SN 2011fu 

The detection of very early time light curve features of 
SN 2011fu, similar to that seen in case of SN 1993J, indi- 
cates a very young age at the time of discovery. Very sharp 
rise followed by a relatively fast decline are explained as the 
detection of cooling phase and depends mainly on the 56 Ni 
mixing and the progenitor radius, as shown by hydrodynam- 
ical models of H-stripped CCSNe (IShigevama et all 1 1994 
Wooslev et all 1 19941 ; iBlinnikov et all ll99Sl ; iBersten et al 



20121 ). For example, in the case of SN 2011dh, for progeni- 
tor radius of < 300 Rq, the cooling phase ends at ~ 5 da ys 
after the explosion (see Figure- 10 of lBersten et alll2012l ). 

In the literature, the first detection of SN 201 lfu 
has been reported to be 2011 September 20.708 (Z. Jin 
and X. Gap, Mt . Nans han, China). However, according to 
ICiabattari et al] |201ll ). this object was not visible on 2011 
August 10 at SN location, putting a stringent limit to the 
explosion date. We collected following pieces of evidences to 
put a constrain on the explosion date of SN 2011fu. 



3.2 Light curve analysis 

In Fig. [21 we plot the calibrated UBVRI light curves of 
SN 2011fu. The LCs span ~ 175 days after the explosion. It 
is clear from Fig.[2]that the photometric observations of this 
supernova started shortly after explosion, showing the early 
declining phase in all bands, which is possibly related to 
the cooling tail after the shock break - out from an extended 
progenitor envelope |Chevalierl. 1 19921; IWaxman et al. 1, 120071 ; 
sonl . 



IChevalier fc FranssonT 120081 ; iNakar fc Saril . 120101) . The LCs 

of SN 201 lfu are strikingly similar to those of SN 1993J, both 
in the initial and the following phases, exhibiting valley- 
like structures followed by rising peaks in all bands. At late 
epochs the LCs are monotonically decreasing in all bands, 
as expected for expanding, cooling ejecta heated by only the 
radioactive decay of 56 Ni and 56 Co. 

Beside SNe 1993J and 2011dh, SN 201 lfu is the third 
known case among lib SNe to date where all the initial 
decline phase, the rise of the broader secon dary peak and 
the final decline have been observed (although lRoming et al.l 
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Table 4. Epochs of the LC valley (t v ) and the secondary peak (t p ) in days after explosion, and their respective apparent magnitudes 
for SN 2011fu and SN 1993J. 



SN 


Band 


LC valley 


Apparent magnitude 


LC peak 


Apparent magnitude 






t v (days) 


clt t y 


t p (days) 






U 


15.51±4.34 


17.67±0.42 


22.93±3.64 


17.43±3.34 




B 


13.75±1.47 


17.93±0.80 


23.29±2.89 


17.51±0.83 


2011fu 


V 


12.87±1.69 


17.45±1.32 


24.96±2.01 


16.95±0.42 




R 


12.95±1.81 


17.01±1.20 


26.40±2.90 


16.50±0.11 




1 


13.50±1.89 


16.86±0.67 


26.64±2.80 


16.41±0.32 




U 


10.33±1.52 


11.94±0.76 








B 


8.82 ±3.36 


12.27±1.74 


19.92±0.70 


11.40±0.17 


1993J 


V 


8.96 ±1.41 


11.89±1.14 


21.67±0.66 


10.87±0.12 




R 


8.81 ±1.06 


11.47±0.61 


22.53±3.24 


10.52±0.37 




I 


9.17 ±1.58 


11.25±0.93 


23.06±1.91 


10.39±0.21 



(2010) reported similar observations for SN 2008ax). In the 
followings, we refer the first minimum of the LC (when the 
initial decline stops and the rise to the secondary maximum 
starts) as the "valley". 

To determine the epochs of the valleys (t v , in days), 
the subsequent peaks (t p , in days) and their corresponding 
brightness values, we fitted a third-order polynomial using a 
X 2 minimization technique to the LCs of both SN 2011fu and 
SN 1993J. The errors in the fitting procedure were estimated 
by the error propagation method. We have taken 1993 Marc h 
27.5 as explosion date for SN 1993J (| Wheeler et all 1 19931 ). 
The derived values of t v , t p and corresponding brightness 
values for both SNe are listed in Table 0] 

The values of t v and t v for both the SNe are similar 
within the errors in all the bands. However, for both SNe, the 
light curves peak earlier in blue bands than in the red bands 
(see Table 0| which is a common feature seen in CCSNe. 
By applying linear regression method, the decline and rising 
rates (in mag day -1 ) were also estimated for three phases, 
i.e. the pre- valley (ai), valley-to-peak (0:2) and after-peak 
phases (013). The results of the fitting are shown in Table [5] 
These values suggest that for SN 2011fu the pre-valley de- 
cay rates (cvi) are steeper (i.e. the decay is faster) at shorter 
wavelengths. This is also true for SN 1993J, where the decay 
rates (ai) were even steeper. Thus, the initial LC decay of 
SN 1993 J was steep er than that of SN 201 lfu during this 
early phase (see also lBarbon et alTll995l ). Between valley to 
peak phase (0:2), the LC of SN 2011fu evolved with a similar 
rate in all the bands, but slower than that seen in case of 
SN 1993J. During the post-peak phase, the LCs gradually 
became flatter at longer wavelengths (see Q3 values in Ta- 
ble [S]). This trend has also been observed for SN 1993 J and 
other Type lib SNe. The B-band LC of SN 201 lfu between 
50 and 100 days after explosion mig ht even show a pla teau, 
similar to SNe 1993 J (see Fig. 3 of iLewis et all 1 1994) an d 
1996cb (see Fig. 2 and the discussions of Qiu et all " 19991 ). 
The plateau-like behaviour of the [/-band LC of SN 201 lfu 
event is more prominent than the (7-band LC of SN 1993J. 

We also determined the Amis parameter for the V- 
band LCs of both SNe, Amis is defined as the decline in 
magnitude after 15 days post-maximum. We got Amis(V) 
= 0.75 mag for SN 201 lfu which is slightly lower than that 
for SN 1993J (Ami 5 (V) = 0.9 mag). Both of these values 




40 60 80 100 
Days since explosion 

Figure 3. Colour curves of SN 2011fu and other Type lib SNe. 
Bottom panel: B — V colour evolution for SNe 2011fu, 2011dh, 
2008ax, 1996cb (symbols) and 1993 J (blue line). Middle panel: 
V - R colour for SN 2011fu and SN 1993J. Top panel: the same 
as below but for the V — I colour. 



are consistent w ith the mean Ami sfy) 
Type Ib/c SNe (|Drout et all l201ll ). 



0.8 ±0.1 mag for 



3.3 Colour evolution and reddening towards 
SN 2011fu 

In Fig. O we compare the evolution of the optical colour 
indices of SN 2011fu with those of other Type lib SNe. 
While constructing the colour curves, we interpolated the 
measured data points (listed in Tabled wherever necessary. 
Before plotting the colours, reddening corrections were ap- 
plied to all the bands. E(B — V) = 0.068 mag was adopted as 
the reddening due to Milky Wa y interstellar matter (ISM) 
in the direction of SN 201 lf u Ischlegel et all 1 1 19981). The 
empirical correlation given bv lMunari fc Zwitterl (| 19971 ) was 
used to estimate the SN host galaxy extinction based on the 
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Table 5. Magnitude decay rate (in mag day 1 ) before valley (cti), rising rate between valley to peak (02) and decay rate after the peak 
{a 3 ) in SN 2011fu and SN 1993J. 



SN 


Band 


Decay rate 


Rising rate between 


Decay rate 






before valley 


valley to peak 


after peak 






K) 


(a 2 ) 


(«3) 




U 


0.24 ± 0.05 


-0.04 ± 0.01 


0.13 ±0.01 




B 


0.15 ± 0.02 


-0.05 ± 0.01 


0.10 ±0.01 


2011fu 


V 


0.11 ± 0.03 


-0.06 ± 0.01 


0.05 ±0.01 




R 


0.09 ± 0.02 


-0.05 ± 0.01 


0.04 ±0.01 




1 


0.09 ± 0.02 


-0.03 ± 0.01 


0.02 ±0.01 




U 


0.38 ± 0.03 








B 


0.24 ± 0.02 


-0.08 ± 0.01 


0.11 ±0.01 


1993J 


V 


0.24 ± 0.01 


-0.10 ± 0.01 


0.06 ±0.01 




R 


0.20 ± 0.01 








1 


0.16 ± 0.01 


-0.09 ± 0.01 


0.05 ±0.01 



measured Nal D lines. For this purpose we calculated the 
weighted equivalent width (EW) of the un- resolved Nal D 
absorption feature in the three spectra (taken on 2011 Oct 
01, 14 and 31, see the log in Table 3), resulted in EW (Nal D) 
- 0.35 ±0.29 A. This corresponds to E(B -V) ~ 0.15 ±0.11 
mag a ccording to the relation given by iMunari fc Zwitterl 
l|l997l ). Finally, we adopted the sum of the two components, 
total E(B — V) = 0.22 ±0.11 mag as the reddening in the 
direction of SN 2011fu. 

The bottom panel of Fig. [3] shows the B — V colour 
evolution of SN 2 011fu along with tha t ofSNe 1993J 
Lewis et all 1 19941). 1996cb dQiu et all 1 19991) . 2008ax 



tudes of SNe Ibc (Type lb, Ic and Ib/c) and II were derived 
using LOSS samples and the average absolute magnitudes 
(close to fl-ba nd as claimed by authors, see discussions of 
iLi et alj (|201ll )) were found to be -16.09 ± 0.23 mag and 
— 16.05 ± 0.15 mag for SNe Ib c and II respectively. In a sim- 
ilar study. [Drout et all |201ll ) also reported that 7?-band ab- 
solute magnitudes of SNe lb and Ic peak arround — 17.9±0.9 
mag and —18.3 ± 0.6 mag respectively. 

Fig. [4l shows the comparison of V^-band absolute 
LC of SN 201 lfu along wit h seven oth er well-observed 
Type lib SNe i.e. 1993J jLewis et all 
' Qiu et all 1 19991 ) 2003 bg ilHamuv et al ' 



Pastorello et all I2OO8I ) and 201 ldh jVinko et alj|20l3 ). It jPastorello et all |200S|). 2009mg (|Oates et all 120121) , 



19941) 
20091) 



1996cb 
2008ax 



is seen in Fig. [3] that the colour curves of SN 201 lfu are 
similar to those of the majority of well-observed Type lib 
SNe, except SN 2011dh which looks being redder than the 
others. 

Similar to SN 1993J, the initial B - V colour of 
SN 2011fu increased (reddened) during the first 10 days 
(note that during the same phase SNe 2008ax and 1996cb 
showed the opposite trend). Between days +10 and +40, 
the B — V colour continued to redden, then after day +40 
it started to decrease and became bluer until the end of the 
our observations. This kind of colour evolution seems to be 
a common trend for Type Ib/c and lib SNe. It may suggest 
that the SN ejecta became optically thin after 40 days. The 
V — R (middle pannel) and V — I (upper pannel) colour 
indices evolve with a similar trend as the B — V colour. 



3.4 Comparison of absolute magnitudes 

The distribution of absolute magnitudes of CCSNe provides 
us infor mation about thei r prog enitors and explosion mech- 
anisms. iRichardson et al.l (|2002l ) made a comparative study 
of the distribution of peak absolute magnitudes in the B- 
band (Ms) for various SNe. They found that for normal and 
bright SNe Ib/c, the mean peak Ms values are —17.61 + 0.74 
and —20.26 ± 0.33 mag, respectively. The Ms values were 
found to be -17.56 ± 0.38 mag and -19.27 ± 0.51 mag for 
normal and bright Type II-L SNe, while for Type II-P and 
Iln SNe the Ms values were found to be —17.0 ± 1.12 mag 

and —19.15 ± 0.92 mag, r espective ly. 

In a recent study by ILi et al.l l|201ll ) , absolute magni- 



2011d h (|Vink6 et all l2012h and 201 lei (|Milisavlievic et all 
|2012| ). For SN 2011fu, the distance D = 77.9 ± 5.5 Mpc has 
been taken from the NEE0 along with total E{B-V) = 0.22 
±0.11 mag as discussed in previous subsection. However, all 
other LCs presented in the figure have been corrected for 
interstellar extinctions and distance values collected from 
literature. Fig. |3J illustrates that the peak My- for various 
Type lib SNe has a range between ~ —16 mag and ~ —18.5 
mag. In this distribution, SN 2011fu is the brightest from 
early to late epochs with a peak absolute magnitude of Mv 
~ -18.5 ±0.24 mag. 



4 BOLOMETRIC LIGHT CURVE 

4.1 Construction of the bolometric light curve 

The quasi-bolometric lightcurve ( UBVRI) was computed by 
integrating the extinction-corrected nuj|f| in all 5-bands. The 
data were interpolated wherever it was necessary and total 
UBVRI flux was integrated using a simple trapezoidal rule. 

In Fig. [5j we compare the quasi-bolometric LC of 
SN 201 lfu along with other three Type lib events i.e. SNe 



7 http://ned.ipac.caltech.edu/ 

8 Fluxes were corrected for interstellar reddening using 
the IDL program ccmjanred.pro available at ASTROLIB 
(http://idlastro.gsfc.nasa.gov/ftp/) by adopting E(B — V) = 0.22 
mag for the total (Milky Way plus in-host) reddening and by as- 
suming reddening law for the diffused interstellar medium {R v = 
3.1). 
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Figure 4. The My light curve of SN 2011fu is compared with 
those of other similar lib events: SN 2011ei, 2011dh, 2009mg, 
2008ax, 2003bg, 1996cb and 1993J. 



Figure 5. The bolomctric light cu rve of SN 2011fu co mpared with 
similar Type lib events SN 1993 J ijLewis et al.Lll994h . SN 2008ax 
iPastorello et all |200sT) and SN 2011dh (Ergon et al. 2012) 



1993J (iLewis et all 1 19941). 20 08ax (jPastorello et all l2008h 
and SN 2011dh (|Ergonl I20121 ). It is obvious that the shape 
of the quasi-bolometric LC of SN 2011fu is similar to that 
of SN 1993J. However, SN 2011fu is more luminous in com- 
parison with the other SNe during the observed phases. 

The un-observed part of the bolometric LC in the Infra- 
red was approximated by assuming blackbody flux distribu- 
tions fitted to the observed R- and /-band fluxes for each 
epoch. At first, we used the Rayleigh- Jeans approximation 
for the fluxes redward of the /-band, and integrated the flux 
distribution between the /-band central wavelength and in- 
finity. This resulted in an analytic estimate for the IR (infra- 
red) contribution as Lir « Xi ■ Fi/3, where Xi and Fi are 
the /-band central wavelength and monochromatic flux, re- 
spectively. Second, we fitted a blackbody to the R- and I- 
band fluxes at each epoch, and numerically integrated the 
fitted blackbody flux distributions from the /-band to radio 
wavelengths (~ 1 mm). These two estimates gave consistent 
results within a few percent, which convinced us that they 
are more-or-less realistic estimates of the IR-contribution. 
Because the /?-band fluxes may also be affected by the pres- 
ence of Ha, we adopted the result of the first, analytic es- 
timate as the final result. The comparison of the integrated 
UBVRI- and IR-fluxes showed that the IR-contribution was 
~ 20 percent at the earliest observed phases, but it increased 
up to ~ 50 percent by day +40 and stayed roughly constant 
after that. 

4.2 Bolometric light curve modelling 

The bolometric light curve (see subsecti on 4.1) was fi t ted by 
the semi- analytic light curve model of lArnett fc Ful (1 19891 ) 
(see also IChatzopoulos et all 120091 ). This model assumes 
homologously expanding spherical ejecta having constant 
opacity, and solves the photon diffusion equation taking into 
account the laws of thermodynamics. This approach was 



first introduced bv lArnettl <ll980T) and lArnettl l| 19821 ). and 
further extended by Arnett fc Ful (|l989T l by taking into ac- 



count the rapid change of the opacity due to recombination. 
The extended diffusion-recombination model was succesfully 
applied to describe the observed light curve of SN 1987A as- 
suming realistic physical parameters |Arnett fc Ful [T989). 

The bolometric LC of SN 2011fu is qualitatively similar 
to that of SN 1987A, because of the presence of the rapid 
initial decline and the secondary bump, after which the LC 
settles down onto the radioactive tail due to the 56 Co-decay. 
This early LC decline in not unusual in Type lib SNe (how- 
ever, see Fig. [2] at early epoch where we compare the LCs 
of SN 201 lfu with SN 1993J), and it is usually modelled 
by a two-component ejecta configuration: a dense compact 
core and a more extended, low er density envelope on top 
of the core <|Bersten et al.Ll2012h . The fast, initial decline is 
thought to be due to the radiation of the cooling outer enve- 
lope (which was initially heated by the shock wave passing 
through it after the explosion), while the secondary bump is 
caused by the photons diffusing slowly out from the inner, 
denser ejecta which is mainly heated from inside by the ra- 
dioactive decay of 56 Ni -» 56 Co -> 56 Fe. After the secondary 
maximum, the decline of the LC is faster than the rate of 
the radioactive decay, which may be due to a recombination 
front moving inward into the ejecta, similar to the condition 
at the end of the plateau phase in Type II-P SNe. 

In order to simulate this kind of LC behavior, we 
sli ghtly modified the or iginal diffusion-recombination model 
of lArnett fc Ful (|l989T l. Instead of having a H-rich, one- 
component ejecta, we added an extended, low-density, pure 
H envelope on to p of a denser, He-rich core. Following 
lArnett fc Ful (|l989l) , we also assumed that the opacity is due 
to only Thompson-scattering, and it is constant in both the 
envelope and the core. Because the envelope was thought 
to contain only H, ft = 0.4 cm 2 g _1 was selected as the 
Thompson-scattering opacity for this layer, while k = 0.24 
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Table 6. Log of parameters derived from bolometric light curve modelling, discussed in section 4.2. 



Parameters 


He- core 


H-envelope 


remarks 




Rprog (cm) 


2 X 10 11 


1 x 10 13 


progenitor radius 




M ej (M ) 


1.1 


0.1 


eject a mass 




k t (cm 2 g 1 ) 


0.24 


0.1 


Thompson scattering 


opacity 


M Ni (M ) 


0.21 




initial nickel mass 




E kln (10 51 erg) 


2.4 


0.25 


ejecta kinetic energy 




E th (0) (10 51 erg) 


1.0 


0.3 


cjccta initial thermal 


energy 
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Figure 6. Comparison of the observed bolometric LC (dots) with 
the best-fit two-component diffusion-recombination model. The 
red and green curves show the contribution from the He-rich core 
and the low-mass H-envelope, respectively, while the black line 
gives the combined LC. 

cm 2 g _1 was applied for the inner region to reflect its higher 
He/H ratio. 

The sy s tem of differential equations given by 
lArnett fc Fvj l|l989l ) were then solved by simple numerical 
integration (assuming a short, At = 1 s timestep which was 
found small enough to get a reasonable and stable solution) . 
Because the photon diffusion timescale is much lower in the 
envelope than in the core, the contribution of the two re- 
gions to the overall LC is well separated: during the first 
few days the radiation from the outer, adiabatically cool- 
ing envelope dominates the LC, while after that only the 
photons diffusing out from the centrally heated inner core 
contribute. Thus, the sum of these two processes determines 
the final shape of the LC. 

Because of the relatively large number of free parame- 
ters, we have not attempted a formal x 2 minimization while 
fitting the model to the observations. Instead, we searched 
for a qualitative agreement between the computed and ob- 
served bolometric LCs. The parameters of our final, best- 
fit-by-eye model are collected in Table [6] while the LCs are 
plotted in Fig. [5] 

It is seen that the best-fit model consists of a dense, 1 
Mq He-rich core and a more extended, low-mass (0.1 M@) 



H-envelope. Thi s is very simila r to th e progenitor configu- 
ration found by iBersten et al.l (|2012l ) when modelling the 
LC of another Type lib event, SN 2011dh, although they 
assumed a more m assive (~ 3 Mp) He- core. Nevertheless, 
it was concluded bv lBersten et al.l l|2012l ) and confirmed by 
the present study that the secondary bump is entirely due 
to radiation coming from the dense inner core of the ejecta, 
and the outer extended envelope is only responsible for the 
initial fast decline of the LC. The estimated ejecta mass for 
SN 2011fu, ~ 1.1 Mq is consistent with the observed rise 
time (~ 24 days) to the secondary maximum of the LC (see 
Eq.10 of IChatzopoulos et all [20121 ). The parameters in Ta- 
ble [6] are also qua litatively similar to the ones derived by 
lYoung et all (|l995l ) for modelling the LC of SN 1993J. 

There are a number of caveats in the simple diffusion- 
recombination model used above, which naturally limit the 
accuracy of the derived physical parameters. The most ob- 
vious limitation is the assumption of constant opacity in 
the ejecta. The pre-selected density and temperature pro- 
files in the ejecta (assumed as exponential functions) are 
also strong simplifications, but they enable the approximate, 
semi-analytic treatment of the complex proble m of radia- 
tive diffusion, as shown bv lArnett fc Ful (|l989l ). Thus, the 
parameters in Table [6] can be considered only as order-of- 
magnitude estimates, which could be signific antly improved 
by m ore sophisticated modelling codes (e.g. IBersten et all 
120121) . 



5 SPECTRAL ANALYSIS 

Properties of the SN 2011fu ejecta were investigated 
with t he multi-pa r ametr ic resonan ce scattering code 
SYH0W jFisher et all Il997l) (see also IBranch et all |2002| ; 
iBaron et all 120051 ; lElmhamdi et all 120061 1 - The~evolution 
of temperature and velocities of layers were traced through 
several months of spectral observations. The SYN0W code is 
based on several assumptions: spherical symmetry; homol- 
ogous expansion of layers (v ~ r); sharp photosphere pro- 
ducing a blackbody spectrum and associated with a shock 
wave at early stages. 

5.1 Comparision between observed and synthetic 
spectra 

In the photospheric phase the spectral lines with P Cygni 
profiles are formed by resonance scattering in a shell above 
the optic ally thick photospher e which produces the contin- 
uum (see IBranch et ail . 1200 ll) . On the other hand, during 
the nebular phase the ejecta is transparent (optically thin) 
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Figure 7. The evolution of the SN 2011fu spectra (grey thick curves, smoothed by a 20A-wide window function) overplotted with SYNQW 
models. The main models are shown by the solid black line. The models with Hj3 fitting are shown as dashed black line. The most 
conspicuous ions are marked. Atmospheric lines are marked with "+". 



in the optical wavelength range. In this case the spectrum 
is dominated by strong emission features including forbid- 
den lines. Each of these two phases of SN evolution can 
be explained with individual approximations and the mod- 
elling of the observed spectra should be made with different 
synthetic codes. There is no sharp boundary between these 
two phases. No strong transition to the nebular phase with 
conspicuous emission features can be seen in the observed 
spectra of SN 2011fu (Fig. [7}. The shape of lines remains the 
P Cyg profile, which suggests that they are formed by res- 
onance scattering, as assumed in SYN0W. Thus, we modelled 
all spectra of SN 2011fu with this code. Before modelling, 
all spectra have been corrected for redshift (see Sect. [5} • 

The strong emission component of the Ha line (prob- 
ably with the C II and Si II contamination) can not be fully 
fitted in the terms of the SYN0W code. We focused primarily 
on the absorption parts of P Cyg line profiles which provide 
information about the expansion velocities of different line- 



forming layers. The SYN0W code allows the usage of different 
optical depth (i.e. density) profiles. Two of them are the 
exponential profile with the parameter of e-folding velocity 
u v e " (r oc exp(—v/v e )) which can be adjusted for each ion, 
and the power-law profile with the index "n" (r oc v~ n ) 
which is applied to all ions in the model. We checked both 
cases and found that the exponential law is more suitable 
for our spectra. The original paper of the SYNQW develop- 
ers and further studies showed a possibility of spectral fea- 
tures which can be detached or undetached from the pho- 
tosphere. These two configurations produce different shapes 
for the line profiles , which was described in the paper by 
ISonbas et alt (|2008h . 

The first three observed spectra are separated by only 
one and two days. That is why they can be modelled by 
similar sets of parameters (see Table 0. Even the spec- 
trum obtained on Oct 14 has a similar continuum slope 
(Tbb ~ 6500 — 6700AT). To verify the pseudo-photospheric 
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temperature derived by SYNOW modelling, we also evaluated 
th e colour temperature (T co i) o f the SN using the model s 
of iDessart fe Hillierl (|2005l > and iBersten fe Hamuvl (|2009h . 
We used the B — V colours for those epochs where spectra 
were available, and then estimated the temperature from 
the corresponding B — V colour. Both of these temperature 
estimates seem to be consistent except for the spectra taken 
on Oct 14 and Dec 22, 2011. 



5.2 Velocity of the pseudo-photosphere 

The velocity of the pseudo-photosphere (an optically thick 
layer, the surface of last scattering for continuum photons) 
can be located by velocities of heavy elements such as Fe II 
and Ti II, which may produce optically thin spectral features. 
However, during the very early phases these features are very 
weak and blended. Therefore, fitting the first three spectra 
by these ions gives a wide range of possible photospheric ve- 
locities, extending from 13 000 to 19 000 km s _1 . The most 
prominent, narrow absorption feature in these spectra is the 
feature near 5650A produced by He I (which may be blended 
with Nal D). This feature is useful to better constrain the 
velocity at the pseudo-photosphere, and decrease the uncer- 
tainty of this parameter at the earliest phases. All velocities 
derived this way are shown in the Vphot column of Table [7] 

5.3 Hydrogen and the 6200A absorption feature 

The wide absorption feature near 6200A can be fitted with 
the help of a high-velocity H-layer (up to V ~20 000 km s _1 ) 
which may be detached from the pseudo-photosphere. On 
the other hand, fitting the emission peak of Ha with SYNOW 
needs lower velocities, but those models cannot reproduce 
the absorption profile (see Fig. [7]). In the V(Ha) column of 
Table [Jj we list the results from the latter, more conservative 
solution. 

The broad absorption at 6200A can be also explained 
with the presence of the C n ion having a high- velocity al- 
most identical to that of Ha. Moreover, Cn also produces 
a small feature near 4400 A. This feature can constrain the 
reference optical depth (r) for ionized carbon. But the con- 
tamination from heavy elements in the blue region makes 
the fitting of the C II 4400A feature uncertain. Thus, the 
presence of carbon cannot be confirmed from these spectra. 

It is also possible to explain the 6200A feature by singly 
ionized silicon. In this case the velocity of Sill must be very 
low. On the other hand, it is expected that the velocity of 
Sin should be equal or only slightly higher than the pho- 
tospheric velocity. It turned out that only the blue wing of 
this wide feature can be fitted by Sill. Although the small 
absorption near 5880 A might be explained by the presence 
of Sin, the observed shape of the 6200A feature does not 
confirm this hypothesis. 

In order to look for other possibilities, we also checked 
different blends of H, Si II, C n and some other ions with dif- 
ferent velocities (assuming undetached as well as detached 
line formation) in our models. At the early phases the range 
of derived velocities turned out to be wide due to the lack 
of observable spectral features formed close to the photo- 
sphere, as discussed above. At the late phases, the wide 
absorption near 6200A splitted into at least three separate 



Table 7. Velocities of the pseudo-photosphere, Ha and H(3 for 
different epochs of SN 201 lfu evolution, derived with SYNOW. 
We assumed that the photospheric velocity {Vphot} ls equal 
to the velocity of Fen. All velocities are given in km s _1 . 
T a b is the blackbody temperature of the pseudo-photosphere in 
Kelvin degrees. The colour temperature (T eo ; ) derived from effec- 
tive temperature — colo ur relations (see lDessart fe Hillierl [20051 ; 
IBersten fe Hamuvl , 120091) is given in last column. 



UT Date 

(yyyy/mm/dd) 


Vphot 

(Fell) 


V(Ha) 


V(H/3) 


T b 


T co i 


2011/09/28 


14000 


16000 


16000 


6700 


6952 


2011/09/29 


14000 


16000 


14000 


6500 


6476 


2011/10/01 


11000 


15000 


11000 


6500 


6052 


2011/10/14 


8000 


11000 


8500 


6700 


5791 


2011/10/29 


6200 


9500 


9500 


5000 


4698 


2011/10/31 


6000 


9000 


9000 


5000 


4795 


2011/11/23 


6000 


9000 


9000 


5000 


5094 


2011/12/22 


5000 


9000 


9000 


5000 


5718 



features (6100A, 6200A, 6350A). These features might be 
explained as a line formation effect for Ha in layers with 
different velocities or the appearance of blending due to ions 
mentioned above. Unfortunately, no firm conclusion can be 
drawn based on the simple parametric models that SYNOW 
can produce. 

The deep absorption near 4700A can be naturally ex- 
plained identifying it as the Hf3 line. We fitted this line and 
the Ha line independently because they cannot be mod- 
eled bv_t4«!_sjmie setof parameters: r, v p hot and v e (see 
also lQuimbv et all [2007). From Table [7J it is visible that for 
the spectra obtained before Oct 29 the fitting of Ha needs 
higher velocities than the fitting of H/3. Although both the 
Ha and H/3 velocities declined in time, the formation of 
the absorption component of Ha stayed at higher velocities 
than for Hft. This may suggest that Ha remained optically 
thick for a longer time than H/3 in the expanding, diluting 
H-rich envelope. 

5.4 Other ions 

To fit the main features in the spectra, we included the 
following elements and ions in SYNOW: H, Fell, Tin, Nal, 
Sill, Ol and Can. Some models were also computed con- 
taining the following elements and ions as alternatives: C n, 
He 1, Fei, Til, Sen, Mgi and the consis t ency of these ions 
were cross-examined with lHatano et all l[l999l '). The heav- 
ier atoms/ions should have velocities close to v p hot but the 
lighter ones may be detached due to e.g. stratification of 
elements in the ejecta. 

In the followings, we show some possibilities to explain 
these features in the observed spectra. Weak features near 
4810A and 6370A can be explained by the presence of low- 
velocity He 1. He can also be found as a blend with Nal in 
the deep absorption near 5650A mentioned above, and as a 
blend with Tin and Fen near 4300A. However, the velocity 
of He 1 may be higher than the photospheric velocity. Even 
in this case, the presence of He can explain all these features. 

A small absorption in the blue wing of Ha at 6630A 
could be modeled by He I or low- velocity Cn. The feature 
near 5050 A could be fitted with Mgi as well. The Can H+K 
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Figure 8. The evolution H a, Hj3 and Fe II line velocities by fitting SYMDW model (see Table [3 ■ Photospheric velocities for SN 2011fu, 
2003bg llHamuv et al.l,l2009h , 1993 J llLewis et al.Lll994lBa.rbon et al.l , ll995l) and 2008ax JPastorello et aUl2008h are shown. The symbols 
of SN 2011fu are connected with lines, those of other SNe with dotted lines. 



feature cannot be fitted well around 3730A, because this 
regime is at the blue end of our observed spectra and all 
of them are very noisy at these wavelengths. But the ab- 
sorption feature near 8400A is compatible with the Can 
IR-triplet. 



5.5 Results of spectral modelling 

Almost all spectral features are described well with elements 
and ions which are usually applied in the case of Type lib 
SNe. However, the strong emission of Ha dominating during 
intermediate epochs can not be fully fitted with the models 
applied. 

The fitting of redder and especially bluer parts has some 
uncertainties due to strong blending of metal lines such as 
Call, Till, Fell and others. Even without precise modelling, 
all spectral sequences can be divided into two groups: the 
first four spectra (up to Oct 14) which are fitted with models 
with the blackbody temperature Tbb ~ 6500 - 6700 K and 
the following four spectra with Tbb ~ 5000 K. 

Generally, the modelling of the SN 2011fu spectra shows 
the decline of the photospheric velocities up to ~ 40 days 
after the explosion (see Fig. [8j). Then, all velocities remain 
approximately at the same, stable level. This behaviour was 
also described in previous works on CCSNe dBranch et all 
2002i: lQuimbv" et al 1. 120071 : iMoskvitin et al ]. l2010h . InFiglgl 
we plot the velocities of Ha, R/3 and Fell for SN 2011fu, SN 
2008ax, SN 2003bg and SN 1993J, illustrating this effect. 



6 METALLICITY-BRIGHTNESS 

COMPARISION OF HOST GALAXIES 

In several earlier studies of CCSNe hosts, it has been already 
mentioned that various SNe subtyp e s occur in different en- 
vironments ( see Prieto et all 2008| : Anderson et all 201C : 
Arcavi et all l2010l: iModiaz et all l201ll : iKellv fe Kirshnerl. 



20121 : ISanders et all 120121 ). Metal licity is a ke y factor in 



all these stu d ies. R ecent studies by I Arcavi et alJ (|2010l ) and 
IPrieto et all (|2008T ) found that SN I b/c host galax i es are 
metal-rich as compared to SN II hosts. iModi az et all (|201ll ) 
found that SNe Ic are more metal-rich (up to 0.20 dex) 
than SNe lb. I n a similar study on SNe Ib/c locations 
iLeloudas et alJ (|201ll ) found a smaller gap between the two 
metallicities (the environment of SNe Ic is richer by ~ 0.08 
dex than lb). In a recent study with a different approach (us- 
ing local emission-line for metallicity esti mates), where 74 
Hn re gions in CCSNe hosts were analyzed, lAnderson et al.l 
i|2010h did not find difference between the metallicities of 
these two environments. 

Type lib host galaxies have been claime d to be more 
meta l -poor than those of SN e lb or Ic (see lArcavi et all 
l2010l : IKellv fe Kirshnerl l2012h . although in another recent 
study which is based on the SN sample from untargeted 
searches (although w ith a rather small sample of 8 SNe lib) 
ISanders et alJ (|2012l ) found that the median metallicity of 
both SNe lb and lib host galaxies is very similar. 

In an attempt to understand the metallicity sce- 
nario for the SN 2011fu host galaxy, we collected the 
latest sample of metallicity data for hosts of CCSNe, 
and their absolute magnitudes in the £?-band from litera- 
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Tremonti et 01. 2004 
Richer & mcColl 1995 
Skillmon et al. 1989 
Kobulnicky & Zorltsky 1999 
Contlni et ol. 2002 
Melbourne & Salzer 2002 



Figure 9. Mctallicity-luminosity relation for various types of 
SN e host galax i es. Th e tiny dots belong to all galaxies used 
by iPrieto et al.l (1200811 (This c atalog is based on SDSS DR4 
lAdelman-McCarthv et al.Ll200r3 . database). Red squares are Type 
II, stars are Type Ib/c and black dots are Type lib SNe, respec- 
tively. The analytic relations collected from several papers (see 
text) are also over-plotted. SN 201 lfu host metallicity is denoted 
by a black triangle. 



tures l|Leloudas et al.l . l201ll ; [Modiaz et al] . l201ll ; [stoll et all 



120121 ) and available online 0. In Fig. [9] the data for these 
host galaxies (36 for lb, 15 for lib and 167 for remain- 
ing II SNe) were then overplotted to the sample contain- 
ing all star for ming galaxies from S DSS DR4 (This sample 
was taken from IPrieto et aD . 120081 ). The relations between 
metallicity and Mb for galaxies from s everal papers are 
also over- plotted Iskillman et all Il989l; iRicher fe McCail 



19951; iKobulnickv fe ZaritskvL Il999l~ IContini et al ., 200 



Melbourne fe Salzeij , |2002| ; iTremonti et allbOoU T 

We estimated the metallic ity of the host o f SN 201 lfu 
using the relation given by iGarnettJ l)2002l ) (see their 
equation 6). We considered Mb = -20.62 mag for UGC 
01626 from HyperLeda. The calculated \og(0/H) + 12 for 
UGC 01626 is 8.901° ™- This value is sli ghtly higher than 
log(0/H) + 12 = 8.55 lArcavi et alll2010t ) and log(0/H) + 
12 = 8.44 jSanders et al.l . l2012l ) for SN type II sample. Our 
analysis using most updated sample of absolute magnitudes 
and metallicities of CCSNe host galaxie s also supports the 
results described in I Sanders et al. (|2012l ). However, it is no- 
ticeable that methods used to determine metallicities are 
based on statistical samples, affected by incompleteness of 
the sample and should be used with caution. 



7 CONCLUSIONS 

We present a comprehensive UBVRI photometeric and low- 
resolution spectroscopic monitoring of the Type lib SN 
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2011fu. To date, only a handful of SNe belonging to this 
class have been observed and studied in detail. 

To the best of our knowledge, our photometric and spec- 
troscopic observations described here are the earliest ones 
reported for this event. The early photometric observations 
strongly suggest the presence of the early-time decline of 
the light curve (which is thought to be related to the shock 
break-out phase) as seen in case of SN 1993J. The early-time 
LC decay rate (qi) of this SN is slower than that derived for 
SN 1993J in all the bands. The rising rates between the LC 
valley to peak (02) observed in SN 2011fu is also somewhat 
slower than in SN 1993J. However, the post peak LC decay 
rate (03) are similar in the two events. 

The colour evolutions of SN 201 lfu were studied using 
our UBVRI band observations. Our data showed that dur- 
ing the very early phases the B — V colour was very similar 
to that in SN 1993J. A similar trend has been found in the 
V — R and V — I colours as well. The evolution of these 
three colours after +40 days were also similar to those seen 
in other CCSNe. The V-band absolute magnitudes of a sam- 
ple of 8 Type lib SNe were compared after applying proper 
extinction corrections and taking into account distances col- 
lected from the literature. In this sample, SN 201 lfu seems 
to be the most luminous event. However, the peak IZ-band 
absolute magnitude of SN 2011fu is not an outlier when is 
compared to the peak brightness of CCSNe of other types. 

The quasi-bolometric LC of SN 2011fu was assembled 
using our UBVRI data and accounting for the IR contribu- 
tion as specified in Section 4. The comparison of these data 
with other known Type lib SNe also shows that SN 201 lfu 
is the brightest Type lib SN in the sample. The bolomet- 
ric LC was modeled by applying a semi-analytical model of 
lArnett fe Ful (|l989l ). This model suggests a 1.1 Mq He-rich 
core and an extended, low-mass (~ 0.1 Mq) H-envelope as 
the progenitor of SN 2011fu, similar to that of SN 2011dh. 
However, the prog enitor radius of SN 2011fu (~ 1 x 10 13 cm) 
turn ed out to be smaller t han that of SN 1993J (~4x 10 13 
cm) jWooslev et all 1 19941 ). The ejected nickel mass for SN 
2011fu was ~ 0.21 M , higher than that of SN 1993J (0.07 
- 0.11 Mq). 

The spectra of SN 201 lfu taken at eight epochs were an- 
alyzed using the multi-parameter resonance scattering code 
SYN0W. The derived parameters describe the evolution of the 
velocities related to various atoms/ions and the variation 
of the blackbody temperature of the pseudo-photosphere. 
The photospheric velocities at the early epochs were higher 
than that of other Type lib SNe. The pseudo photospheric 
temperatures were found to be between 6700 K and 5000 
K, decreasing from initial to later phases. The tempera- 
tures from SYN0W were also checked by comparing them with 
colour temperatures calcu l ated from B — V vs T e ff relations 
|Bersten fe Hamuvl 120091 ; iDessart fe Hillierl . |2005| ). These 
different temperature estimates were found being consistent. 
The appearance of the main observed spectral features were 
also successfully modeled with SYN0W by assuming H, He I 
and various metals (mostly Fen, Tin and Call), which are 
typical of CCSNe spectra. The estimated value of the metal- 
licity of the host galaxy of SN 2011fu is 8.901°,;^ similar to 
those for other Type lib SNe. 
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